Transverse Meissner Physics of Planar Superconductors with Columnar Pins 
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The statistical mechanics of thermally excited vortex lines with columnar defects can be mapped 
onto the physics of interacting quantum particles with quenched random disorder in one less di- 
mension. The destruction of the Bose glass phase in Type II superconductors, when the external 
magnetic field is tilted sufficiently far from the column direction, is described by a poorly under- 
stood non-Hermitian quantum phase transition. We present here exact results for this transition in 
(l+l)-dimensions, obtained by mapping the problem in the hard core limit onto one-dimensional 
fermions described by a non-Hermitian tight binding model. Both site randomness and the relatively 
unexplored case of bond randomness are considered. Analysis near the mobility edge and near the 
band center in the latter case is facilitated by a real space renormalization group procedure used 
previously for Hermitian quantum problems with quenched randomness in one dimension. 



I. INTRODUCTION 

The physical properties of vortices in random pinning 
potentials have been the focus of investigation for many 
years. Recent progress has taken place in a particular 
aspect of this problem - vortices trapped inside a slab 
of superconducting material with columnar defects. In 
general, the problem of vortex pinning in a superconduc- 
tor with columnar defects can be recast via the trans- 
fer matrix method in terms of the non-Hermitian quan- 
tum mechanics of bosons with a constant imaginary vec- 
tor potential. 1,2 In the case of a thin superconducting 
slab, the mapping to non-Hermitian quantum mechanics 
leads to a one-dimensional problem. In this case there 
has been recent progress in the following specific prob- 
lems: (1) Vortex hopping in a regular array of colum- 
nar pins with on-site and nearest-neighbor repulsion. 3,4 
(2) Transmission through a weak- link. 3,4 (3) Vortex tun- 
neling with Cauchy-distributed random on-site pinning 
potential (Lloyd model) 5-7 (4) Some progress in under- 
standing vortex dynamics with on-site and hopping en- 
ergies both random was made in Ref. 8-10. 

In this paper we first concentrate on a fifth case, 
namely, random hopping in an otherwise uniform pin- 
ning array. The approach we take here is the strong 
randomness real-space renormalization group (RSRG). 
This approach, pioneered by Ma, Dasgupta, and Hu, 11,12 
and developed further by Fisher, 13,14 is very successful in 
treating one dimensional random spin chains as well as 
other models. 15-17 Use of this technique in the context of 
non-Hermitian transfer matrices, however, appears to be 
new. We note that recently a related problem of a vortex 
lattice pinned by a single columnar defect was considered 
by Radzihovsky in Ref. 18. 

This problem is of interest for several reasons. 
First, although much is known about non-Hermitian 
tight binding models with site randomness in (1+1)- 
dimensions, 2,6-10 the case of pure bond randomness (cor- 
responding to irregularly spaced columnar pins of equal 
strength in a slab) is relatively unexplored. All states are 
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FIG. 1: Schematic phase diagram of a flux liquid with colum- 
nar pins at a fixed external magnetic field Hu parallel to the 
pins as a function of temperature T and tilted external field 
H±. The transitions along the line H±(T) map onto the non- 
Hermitian quantum phase transition problem discussed in this 
paper. Inset shows the behavior of B± near HI, where the 
transverse Meissner effect breaks down and the vortices begin 
to tilt. 



localized for the Hermitian problem of vortices in a thin 
slab without an external tilt field when only site random- 
ness is present. The tilt field h must then exceed a finite 
threshold before delocalized states appear in the center 
of the band. 2 In contrast, there is always one delocal- 
ized state exactly at the band center for the Hermitian 
problem with only random hopping, accompanied by a 
diverging localization length of the localized states on 
either side. 19 

As we show explicitly in Sec. VII, additional delocal- 
ized states then appear immediately in the center of the 
band for any nonzero tilt in the thermodynamic limit. As 
discussed further below, for larger values of h both prob- 
lems do have similar mobility edges, separating localized 
states near the band edges from delocalized states near 
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the band center. 

Perhaps more important, the analytically tractable 
free fermion model discussed here sheds light on the 
poorly understood non-Hcrmitian quantum phase tran- 
sition that describes the physics of vortices as one at- 
tempts to tilt them away from the direction preferred 
by columnar defects. Suppose for simplicity the longi- 
tudinal applied magnetic field (i.e., H\\, the field par- 
allel to the columnar pins) produces a vortex density 
which docs not exceed the density of column pinning 
sites. Then, as shown in Fig. 1, the low temperature 
Bose glass phase, with essentially all flux lines localized 
on columnar defects, is expected to be stable for a small 
additional field H± perpendicular to the column direc- 
tion. The transverse applied field H± is proportional to 
the tilt field h discussed in this paper. Below the zero 
tilt Bose glass transition temperature Tbg, perpendic- 
ular fields less than a critical value HI leave vortices 
untilted and trapped on columnar pins in the thermody- 
namic limit. 1 The equivalent quantum problem involves 
interacting quantum bosons in a disorder potential and 
constant imaginary vector potential proportional to H±. 2 
Provided the transition is not first order, the breakdown 
of this transverse Meissner effect above H]_ can be de- 
scribed by a critical exponent (, according to 

B ± ex (H± - Hif , (1) 

where B± is the transverse flux due to the tilted vortices. 
Heuristic random walk arguments based on the entropy of 
vortices wandering in the presence of thermal fluctuations 
lead to the estimates, 20 

C=|,(d = 3) 
C=h(d = 2) 

in three and two dimensions respectively. However, there 
are reasons to doubt these predictions. First, as will 
be become clear, at least for the exactly soluble (1+1)- 
dimensional model discussed in this paper, the non- 
Hermitian quantum phase transition induced by tilt ef- 
fectively occurs at a finite wave-vector, calling into ques- 
tion simple random walk arguments based on physics at 
k = 0. Second, although the phase diagram shown in Fig. 
1 is well established experimentally, 2 1-23 the one exist- 
ing experimental measure of the exponent £ 24 (based on 
current-voltage characteristics, see below) is inconsistent 
with Eq. (2). Indeed, the experiments on bulk super- 
conductors with columnar pins in Ref. 24 find ( = 1/2, 
which disagrees with the prediction of Eq. (2) in d = 3. 
In fact, we show in this paper that, for the special case 
of free fermions in (l+l)-dimcnsions (corresponding to 
Luttingcr liquid parameter g = 1 in Fig. 4) one obtains 

C = 1 (free fermions, 1+1 dimensions), (3) 

which disagrees with the prediction of Eq. (2) in d = 2. 

In principle it might be possible to check predictions 
like those in Eqs. (2) and (3) by magnetic torque mea- 
surements, which are sensitive to the differences in the 



direction of B and H that are an essential part of the 
transverse Meissner effect. 22 Alternatively, as discussed 
in Refs. 1 and 20, the exponent £ determines a density 
of kink excitations connecting nearby columnar defects 
(see Fig. 2), which in turn control the linear flux flow 
resistivity above H^_. Current- voltage curves ar highly 
nonlinear below HI, and the linear resistivity vanishes. 1 
To understand the linear resistivity above Hj_, consider 
the geometry shown in Fig. 2a, where a current flows per- 
pendicular to the plane defined by the column direction 
and the average field direction defined by a set of tilted 
vortex lines. Imagine first a bulk sample, of dimensions 
L x W x W, where L is the sample length along the 
columns. Suppose this field is inclined at a small angle 
9 = B±/B\\ away from the column direction. Then a 
typical perpendicular distance Ax traversed by a vortex 
across a sample is Ax — 6L. If the density of columnar 
pins is n p , this tilt leads to 

Nl = Ax/n-^ = ^ n y*L (4) 
B \\ 

kinks associated with a single vortex. In clean Type II 
superconductors, at temperatures high enough so that 
residual pinning by point impurities can be neglected, 
these kinks will slide along columns, due to the Lorentz 
force Jl = 4>$J l c caused by the current, where 4>o is the 
flux quantum and c the speed of light. Upon multiplying 
by the total number of flux lines W 2 B\\/4> Ql where, we 
find a gas of kinks with density 

n k = B^n^/fo (5) 

per unit volume. As expected, rife is proportional to B±, 
which leads using Eq. (1) to a flux flow resistivity 

p = Po B ± e/fr) cc (H ± - Hl) C , (6) 

where po is the normal state resistivity, £ (the coher- 
ence length) is the size of the normal vortex cores which 
contribute to the dissipation, and we have neglected co- 
efficients of order unity. A similar calculation can be 
carried out for a two dimensional slab with dimensions 
L x W x d, with d <C L, W and a single sheet of columnar 
pins along L with spacing n" 1 . We also require d < A, 
where A is the London penetration depth. We then find 
that the flux flow resistivity for a current perpendicular 
to the slab (i.e., along the direction d) is identical in form 
to the three dimensional result (6). 

We hope that the calculations in this paper, although 
only valid in (1+1) dimensions at a special temperature 
deep within the Bose glass phase (see Fig. 4), will stim- 
ulate further theory as well as experiments aimed at a 
more complete determination the exponent £, in both 2 
and 3 dimensions. 

A. Statement of the problem 

When a magnetic field is applied to a planar super- 
conductor with parallel columnar defects in the plane, a 
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FIG. 2: (a) Schematic of a single tilted vortex line in a sample of thickness L and inclined at an average angle 9, interacting 
with a row of variable strength columnar pins (represented by solid dashed and dotted vertical lines) used in our estimate of 
the flux flow resistivity in the presence of a current J flowing into the plane of the diagram, (b) Top view of this same situation, 
showing the path of the tilted vortex projected down the z-axis. n p is the density of columnar pins. 



competition occurs between the tendency of the vortices 
to be pinned by the columnar defects, and the portion of 
the in-plane magnetic field which is normal to the colum- 
nar defects. The resulting lock-in of vortex trajectories 
parallel to the columns up to a critical tilt field is the 
transverse Meissner effect discussed above. 

In this work we concentrate on a model of vortex 
hopping in a random array of columnar defects (see 
Fig. 3). As described previously 1,2 , this statistical me- 
chanics problem is mapped onto a quantum mechanical 
boson-hopping problem. When the magnetic field that 
produces the vortices is tilted relative to the columns, our 
model introduces non-Hcrmiticity to the Hamiltonian. 
The Hamiltonian whose exponential gives the transfer 
matrix is: 



W = E (~Wi (blb.+re-' 1 * + bj +1 b 4 e^) 

+ {a - (j,) bjb, + %m (m - 1)) . 



(7) 



Each vortex is represented by a boson with creation and 
annihilation operators b,b^. As discussed above, we re- 
strict our attention to hard core vortices, and set the 
on-site repulsion to infinity (U — ► oo). 

The average number n of vortices per pinning site is 
equivalent to the longitudinal magnetization induced by 
a magnetic field H\\ parallel to the defects. The part 
of the magnetic field H±, tilted relative to the colum- 
nar defects, induces an imaginary vector-potential hi on 
each bond. We expect it to depend linearly on the dis- 
tance between columnar pins. Note, however, that using 
a similarity transformation one can redistribute the hi's 
such that each bond carries the same imaginary vector- 
potential 8,9 . In the following we will therefore eventu- 
ally apply a uniform tilt h for all bonds. In Eq. (7), \i 
is the vortex chemical potential, which is controlled by 
the external magnetic field and the depth of the pinning 



potential. The hopping energies Wi and pinning ener- 
gies Cj are random variables with some relatively well 
behaved distribution. As is shown schematically in Fig. 
3, the Hamiltonian (7) describes the statistical mechanics 
of vortices hopping between columnar defects with an ex- 
ternal magnetic induction tilted relative to the columnar 
pins. 

In general, our goal is to characterize the transverse 
Meissner effect in the superconductor. For that purpose 
we also need to define the transverse magnetic flux in 
terms of the bosons of Eq. (7). Quite intuitively, the 
transverse flux is equivalent to the boson current, and is 
found by differentiating the Hamiltonian with respect to 
the external transverse field: 1 ' 2 

J * = H) lr = HK ( b ^+ ie ~' is - b ll b * e ' l, ) • 

(8) 

We will consider and contrast two cases: first, the case 
of random pinning energies with uniform spacings, and 
then the case of identical pins but random spacings. We 
will show that these two cases are quite different, and 
proceed to characterize the case of random spacings in 
some detail. 

As we show in Sec. II B, the Hamiltonian (7) with 
the hard-core requirement (U — ► oo) is equivalent to a 
fermion model without interactions. In general these in- 
teractions are quantified by the Luttinger parameter 4 



7tT 



9 



VC11C44 



(9) 



where T is the temperature, Cn and C44 are the vor- 
tex compressibility and tilt modulus respectively. In this 
work we concentrate on the line g — 1, and ignore near- 
est neighbor and higher range interactions between the 
hard-core vortices. The general phase diagram for the 
interacting vortex lines with tilt in (1 + l)-dimcnsions is 
shown in Fig. 4. 1,20 
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FIG. 3: The exponential of the Hamiltonian (7) determines 
the transfer matrix describing the statistical mechanics of vor- 
tices (black lines) fluctuating in a planar superconductor with 
parallel columnar defects piercing it at random intervals (gray 
stripes). The distance between the defects, ctj, as well as their 
energy depth (illustrated by their width), Vj, are random vari- 
ables. The irregular spacing leads to strong randomness in the 
vortex hopping, Wj oc e~ laj . The external magnetic induc- 
tion and its direction relative to the columnar pins determine 
the vortex chemical potential, fi, and the effective tilt hj oc a,j . 
In this paper we compare and contrast the case of identical 
pins with random spacings, and the case of random pinning 
energies but uniform spacings. 



B. Summary of results 

In the first part of this paper (Sees. II and III) we re- 
view and derive properties of the random-pinning Lloyd 
model. We concentrate on the minimum tilt required 
to form delocalized eigenstates, i.e., the critical tilt as a 
function of parallel field to destroy the transverse Meiss- 
ner effect, and the resulting transverse flux as a function 
of both tilt and parallel field. The derivation of these 
properties provides us with a baseline for a comparison 
of the random pinning model with the random-hopping 
model. 

In the second part of the paper we concentrate on 
random-hopping. In this model all pinning sites are as- 
sumed to be identical, but with random distances be- 
tween them (since hopping depends exponentially on the 
inter-site distance, the hopping strengths, Wi, will be 
strongly random). For this purpose we employ the real- 
space renormalization group (RSRG) method. Using this 
method we obtain the density of states for the vortex- 
hopping Hamiltonian. We derive this method for the 
vortex-hopping problem in Sec. IV. 

The random-hopping problem has two unique features 
that are connected to each other: the localization length 
of the vortex eigenfunctions diverges near the middle of 
the band (E — 0), and there is also a singularity of the 
DOS at the same place. We use the DOS to derive a rela- 



Vortex 
3/2/ 


Liquid 




y i 


- - - - V 








/ Bose 


Glass \ 




h 



FIG. 4: The phase diagram for interacting vortices, with effec- 
tive Luttinger parameter g proportional to temperature (see 
Eq. 9), tilt h, at a fixed parallel magnetic field H«. For real 
bosons, the vortex liquid phase for g > 3/2 and vanishing tilt 
h — corresponds to a superfluid with off-diagonal long range 
order. Little is known about the Bose-glass depinning tran- 
sition at finite h. In this paper we investigate this transition 
at o=l. 



tionship between the effective vortex chemical potential, 
fi, and the applied parallel field Bu. We then proceed 
to show that any nonzero tilt will produce delocalized 
eigenstates, and derive the critical tilt h c (b) at which the 
transverse Meissncr effect breaks down. We also derive 
the localization length of vortex states near the mobility 
edge. 

By employing a simple spectral formula explored in 
Refs. 6 and 8,9, together with the results of the real-space 
RG and the general arguments in App. A, we derive the 
following properties: the angle of approach of the de- 
localized spectrum in the complex plane (dlmE / dReE) , 
the contribution to the transverse magnetization of a sin- 
gle delocalized vortex state, and the total vortex current, 
or transverse flux, near the breakdown of the transverse 
Meissner effect. 

We emphasize that the results derived here assume 
that the system is in the universal low energy limit. A 
numerical investigation of this limit requires large sys- 
tem sizes that allow the RSRG to reach low energies. 
This is confirmed in Sec. VII where finite-size systems 
of vortices and columnar pins with random pinning and 
random hopping are diagonalized exactly. 



II. EQUIVALENT MODELS 

Hamiltonian (7) in the limit U — > oo is not as famil- 
iar to us as some other equivalent models. In this sec- 
tion we map the boson Hamiltonian (7) to a spin model 
with easy plane anisotropic interactions (XX), and to a 
fermion hopping model. These mappings will be useful 
when applying the RSRG in Sec. IV and exact diagonal- 
ization for finite systems in Sec. VII. 
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A. Mapping to a spin model 

By the simple transformation of the boson operators 
on the n th lattice site, 

t>t = (-i)» . £+ = (-i)« . (§z + i§y), 

b n = • s- = (-i) n ■ (s% - isy) , (io) 

btb„ = l/2 + ^, 

the Hamiltonian in Eq. (7) is transformed to an XX 
fcrromagnet in an external magnetic field: 



n = E (-2^ ((sfSf +1 + sfsf +1 ) cosh(/i) 

+i (sfSf +1 - SfSf +1 ) sinh(i 
+ (c . 5?)). 



In the spin variables, a vortex pinned at site i corre- 
sponds to = +1/2; an empty site i is transformed 
to a site with Sf = —1/2. The transformation in Eq. 
(11) is the special spin-1/2 case of the Holstein-Primakoff 
transformation. 25 In its higher spin version it has factors 

of the form \fl — bjjb„ which for spin-1/2 reduce to ei- 
ther zero or one. 

At half-filling (i.e., when fi = ej = 0) and zero tilt h the 
ground state of Hamiltonian (11) is the random singlet 
phase. In this phase singlets form in a random fashion 
between sites connected by strong Wi. Once a singlet 
forms, the nearest neighbors of the two sites involved 
also interact, but with suppressed strength. By pairing 
up strongly-interacting sites into singlets, we iteratively 
reduce the energy scale of the Hamiltonian, until we ex- 
haust all sites. At this point, singlets connect many near- 
est neighbors, but occasionally singlets connect very far 
away spins, leading to power-law decaying correlations 
(for a review, see Ref. 13). The idea behind the random 
singlet phase and its formation are shown in Fig. 5. A 
thorough discussion of this phase and its implications for 
vortex pinning will be given in Sees. IV and V. 



B. Mapping to fermion hopping problem 

By using a version of the Wigner- Jordan 
transformation, 26 we can also map the hard-core 
boson Hamiltonian, Eq. (7) with U — oo, to a fermionic 
random hopping problem. This mapping is also quite 
straightforward; we start with 



n— 1 + n— 1 + 

b„^ n e^ c „, bt^ n e-^'ct. (12) 

j=-oo 



n-1 

n 

j = -oo 



where the string operator ensures the anti-commutation 
relations of the c„'s. Now we can write the Hamiltonian 



w, w. 



FIG. 5: (a) In a random XX spin chain, described by Eq. (11) 
with h = 0, strong bonds such as W2 localize a spin-singlet. 
Quantum fluctuations induce a coupling between the neigh- 
bors of the singlet. This coupling has the same XX form, 
but a substantially reduced scale w e ff = W1W3/W2 which is 
much smaller than wi, W2, and 103. By repeating the singlet 
formation process, the energy scale of the effective Hamilto- 
nian is reduced, (b) The random singlet state. Singlets form 
in a random fashion, mostly between nearest neighbors, but 
they also connect very far away spins. Long distance singlets 
give rise to average correlations that decay algebraically with 
distance. 



of the vortices as though they are fermions: 

(13) 

Similarly, the local current operator becomes: 

J« = (^*)|^7 
= (-i)wi (h\b l+1 e- h * - bj +1 b l e /l -) (14) 

= {-i)wi (c\c l+1 er hi - cj +1 c l e /li ) . 

In the absence of the tilt field, the ground state of the 
Hamiltonian (13) is well understood. 13 When the pins are 
identical (e^ = 0), the random hopping localizes all states 
except at half filling, where there is always a delocalized 
state. The localized states away from the mobility edge 
are related to the random singlets in the XX spin chain 
of Eq. (11). Instead of a singlet, however, pairs of sites 
share a single fermion. 

In this model, half filling is obtained when fi — > CP. 
In this limit, the last fermion inserted in the system is 
delocalized between two sites with a distance that is of 
the order of the system size, i.e., it is delocalized. All 
other fermionic states, however, are localized. 

In the presence of the tilt field, h, an entire band of 
delocalized states appears, and the lower mobility edge 
moves down to fillings below one half, and to negative 
chemical potentials fj, < 0. By using the mapping to 
frcc-fcrmions (Eq. 13) and the known real space renor- 
malization group (RSRG) results for this model, we can 
obtain much insight into the delocalized phase. In partic- 
ular, we will derive a universal relationship between the 
field h, the mobility edge fih, and perhaps most impor- 
tantly, the vortex density (i.e. parallel magnetic field) at 
the mobility edge, ph. 
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III. RANDOM PINNING ENERGY AND 
UNIFORM HOPPING - EXACT RESULTS FROM 
THE LLOYD MODEL 

In this section we will review exact results for the 
Lloyd model. 3 ' 5,6 We consider the one-dimensional quan- 
tum Hamiltonian, Eq. (7) with infinite on-site repulsion 

(U->oo): 

n = T,{- w (b+bi+ie- 71 + h\ +1 h ie h ) - ( Ci - /i) bjb, 

(15) 

where bj and are creation and annihilation operators 
of hard-core bosons, which represent the vortices. The 
hopping matrix element w is now site independent. We 
could equally well use the fermion representation in Eq. 
(13) in what follows. We set the lattice constant to 1, and 
consider a lattice of length L sites. For small tilts, the 
imaginary gauge- field h is proportional to the angle of the 
applied magnetic field relative to the columnar defects. 
The parallel applied field translates into the chemical po- 
tential for the vortices via the constraint, 



/ 



g{e)de = n, 



(16) 



where n is the average number of vortices per columnar 
defect, and g(e) is the density of states associated with 
the Hamiltonian of Eq. (15). Assuming that the dis- 
tribution of pinning energies, e i; is symmetric, we can 
invoke particle hole symmetry, and rewrite Eq. (16) as: 

n 

Jg(e)de= J g(e)de - J g(e)de = 0.5 - n. (17) 



In this section we will review the known qualitative 
feature of this non-interacting model of vortices with 
random pinning strength e,. We will also focus on the 
analytic results known for the Lloyd model, 5-7 which is 
described by the Hamiltonian (15) with a special distri- 
bution of the pinning energies 
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This model has played an important role in our under- 
standing of one-dimensional localization of electrons. 



Structure of the spectrum and critical 
derealization tilt 



Much of the physics of interacting vortices in (1 + 1)- 
dimensions can be inferred from the shape of the spec- 
trum of the Hamiltonian as a function of the applied 
transverse field h. For zero tilt h, the entire spectrum of 
Hamiltonian (15) is due to localized states, and is there- 
fore real. 2 The density of states changes as a function 



of energy but experiences no singularities. When the 
external magnetic field is tilted relative to the colum- 
nar defects, the spectrum of the localized states does 
not change, but the wave functions associated with them 
spread out in the direction of the tilt. The invariance of 
the spectrum to the tilt, as long as the wave-functions are 
localized, can be easily understood, since for any eigen- 
state of the Hamiltonian we can carry out a gauge trans- 
formation such that all the imaginary gauge-field in a 
sufficiently large system is on a bond between sites where 
the eigenstatc has no support. 

Above some critical tilt ho > 0, a subset of the eigen- 
functions delocalize, and their energies become complex 
(see Fig. 6). As shown in Ref. 8,9, the parts of the spec- 
trum that first become delocalized are at the maximum 
of the DOS g(e), which is generally in the center of the 
spectrum, at e = 0. 

The transverse Meissner effect, however, persists until 
the vortex-eigenstates at the "Fermi energy" fi become 
delocalized. This Fermi energy for the vortices is deter- 
mined by the parallel (longitudinal) magnetic field By, 
or equivalently by the vortex density n per pin. The tilt 
at which the transverse Meissner effect breaks down is 
thus a function of n, and it is always bigger than the 
threshold-tilt h : 



h c (n) > h . 



(19) 



The qualitative description above is exemplified by the 
spectrum of the the Lloyd model. 5 ' 6 The delocalized por- 
tion of the spectrum of the Lloyd model is known to be: 5 



E^ = — 2wcos(±fc„ + ih) =F ij, 



(20) 



where 7 is controls the width of the site-randomness dis- 
tribution in Eq. (18). As the tilt increases such that 
h > ho, a bubble of delocalized states appears in the 
center of the band (c.f. Fig. 6). When h > h c (n), this 
bubble engulfs the chemical potential. 

From the spectrum of the Lloyd model, Eq. (20), 
one can easily derive ho and h c . The lowest critical tilt 
ho that produces delocalized eigenfunctions is obtained 
when the imaginary part of the energy becomes nonzero, 
lmE n 7^ 0, for some eigenstatc E n . As is obvious from 
Fig. 6b, this derealization tilt will be the critical tilt 
where the transverse Meissner effect disappears at half 
filling (n = 0.5). h is given by: 



sinh(/i ) = 



J_ 

2w 



(21) 



The range of wave- vectors k n that defines the delocalized 
states are thus given by: 



sin IfcJ > 



7 



2w sinh(/i) 



(22) 



To obtain the critical tilt as a function of longitudinal 
magnetic field, we observe that the number of delocalized 
single- vortex states is 



2L ( 7 

Ndei = — arccos . 

tt \2wsmh(h) 



(23) 
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g(E) h=h < 




ImE 



ImE 



ReE 



h>h, 



ReE 



h>h c (n) 



ReE 



(note that the argument of the arccos is never bigger 
than 1). Thus the critical longitudinal vortex density for 
a given tilt is: 



nc = ^^ arCC ° S ( 2w si 7 nh W )' (24) 



which can be inverted to give the critical tilt for a given 
magnetization: 



sinh(/i c 
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B. 



2w cos(7r(0.5 — n)) 
Transverse magnetic flux near critical tilt 



(25) 



Once vortices form delocalized states, the magnetiza- 
tion in the superconductor is no longer parallel to the 
columnar defects, and transverse flux appears. The to- 
tal transverse flux, or vortex current in the quantum- 
mechanical picture, can also be found out easily for the 
Lloyd model. 5 

In App. A we derive a general rule for the total vortex 
current for delocalized states (Eq. A4): 



'total 



dh 



1 



{ImE,, - ImE, 



(26) 



FIG. 6: (a) All eigenfunctions of the random-pinning Hamilto- 
nian are localized and independent of the tilt provided h < ho . 
We suppose that the density of states (DOS) is peaked around 
E — but is not singular. The effective chemical potential 
of the vortices, /i, is determined by the longitudinal mag- 
netization n, such that the density of occupied states below 
fi equals n. (b) When the tilt of the applied field becomes 
h > ho, a bubble of delocalized states develops where the 
DOS is peaked. As long as the bubble does not reach /j, there 
is still no transverse magnetization. The thick lines mark 
the support of the eigen-energies in the complex E-plane. (c) 
The transverse Meissner effect breaks down when h = h c , and 
the delocalized bubble reaches fi. The total vortex current is 
shown in App. A to be proportional to the imaginary part of 
the energy eigenvalue at the Fermi surface (c.f. Eq. A4). 



where ImE^ is the imaginary part of the energy eigen- 
value with real part /i and fj. m (h) is the chemical poten- 
tial above which vortex states are delocalized (mobility 
edge). L is the total number of lattice sites. As discussed 
in Ref. 2,4, this total 'current' (similar to the derivative 
with respect to vector potential which gives the current in 
a quantum system) is proportional to the perpendicular 
component of the magnetic flux. If we consider our sys- 
tem having longitudinal vortex density n and tilt h > h Cl 
then the number of delocalized vortices (as opposed to 
vortex states) is n — n c where n c is given in Eq. (24). 
Therefore, the k n that are occupied by delocalized vor- 
tices are given by 



arcsin n ■ r /,n = fc c < | fc n | < arcsin . + tt (n - n c ) = k F (27) 

\2wsmh{h) J \2w sinh(/i) / 

Now, using Eq. (A4) (see Appendix A) we can directly write the total current as: 

Jtotai = 7 (JmE kp - lmE ka ) = ±2™ sinh(/i) • (sin (arcsin 2ws ? nh(h) +n(n- n e j) - sin (arcsin 2u)Bi ^ h(h) ) ) 

V V ' V ' J (28) 

= \2w sinh(ft) • ( 2wsi ; h(/t) (cos (tt (n - n c )) - I) + ^1 - (2sin g fe)] ^ sin (tt (n - n c ))) . 
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Near the derealization transition we have: 



Jtotai = 2wsmh(h)\ 1 



7 



^2wsinh(/i) 
Using the derivative of Eq. (24) we obtain 
dn c 1 7 cosh(/i) 1 



(n-n c ). (29) 



d/i 27T 2w sinh /i 



^2u> sinh(h) ^ 



Upon putting this back in Eq. (29) with h — h c s=s 
\dn/dh ■ (n — n c ) | we obtain a simple expression for the 
total current in the Lloyd model near derealization, 
namely 



Ji 



total 



2tt tanh(/i) 



(ft- 



(30) 



This linear onset of the vortex current is qualitatively 
similar to numerical results for a box distribution of pin- 
ning energies, as presented in Section VII. When reinte- 
grated as a formula for the transverse magnetic flux of a 
vortex system, Eq. (30) leads to the prediction £ = 1 for 
the exponent discussed in the introduction, Eq. (2). 



IV. REAL SPACE RG OF THE BOSON 
HOPPING PROBLEM 

One purpose of this paper is to understand the prop- 
erties of vortices hopping in an array of identical pins 
with random locations, and contrast these properties 
with those of the random-pinning energy system with 
uniform hopping. The random-hopping problem is de- 
scribed by the Hamiltonian in Eq. (7) with U — > oo, 
ej = 0, and u>i random. In order to solve this model 
we apply the real-space renormalization group (RSRG) 
procedure. 11-14 In this section we review the application 
of the RSRG to the random hopping boson Hamiltonian. 

The RSRG allows us to diagonalize the random Hamil- 
tonian iteratively by eliminating the high energy degrees 
of freedom. The first step is finding the strongest bond in 
the chain and diagonalizing it while ignoring the rest of 
the Hamiltonian. Quantum fluctuations give rise to new 
interactions that bridge over the bond we diagonalize. If 
we are lucky, these renormalized interactions will be of 
the same form as the original Hamiltonian. Let us carry 
this out for the boson- hopping Hamiltonian Eq. (7). In 
effect, we diagonalize the Hamiltonian (7) iteratively by 
putting vortices into a hierarchy of localized states shared 
by pairs of sites. 



A. Diagonalization of the strongest bond 

The Hamiltonian we are interested in is given in Eq. 
7, with U — > oo and = 0: 



Suppose Wi is the strongest hopping energy in the chain. 
We first diagonalize a single nearest neighbor bond: 

Ht = - Wi (h\b l+1 e- h + bj +1 b J e' l )-Mbjb J + r4 +1 b m ) 

(31) 

Because the vortices are hard core in this limit the only 
occupied states available in the Hilbert space are 



|i> = bj| 



|i + l) = bj +1 |0). 



(32) 



Diagonalizing Eq. (31) in the subspace spanned by (32) 
means diagonalizing the matrix: 



— jU — Wie 
h 



(33) 



-Wie 



Note that the eigenvalues of this Hamiltonian can not 
depend on the tilt, in accordance with what we know 
on localized states - they carry no transverse magnetic 
field. 2 The eigenvalues and right eigenvectors of (33) are 



E± = —fj, ± W r 



\±) = -L(e- h / 2 \i) T e h / 2 \i + l)). 



(34) 



Note that in this non-Hermitian quantum problem, the 
left eigenvectors are not the Hcrmitian conjugates of the 
right ones. In fact: 

(±| = -L(e^| Te -^ + l|). 

As long as E_ < and E + > the ground state of Hi 
from Eq. (31) will contain a single vortex split between 
sites i and i + The cigenstate of bond i is therefore 



|_) = (e-V2 b t +e V 2b t +i )| 0) . 



(35) 



This state is described in Fig. 7a. 

The next step is to include the rest of the Hamilto- 
nian, and in particular neighboring bonds corresponding 
to Hi±i. These terms will lead to quantum fluctuations 
above the ground state of Hi. Thus it is necessary to 
find the allowed excited states of (7) which involve 2 and 
vortices. These states are simply 



IM) 



btb| +1 |0>, 



|0,0> = |0>, 



(36) 



with energies E = —2{i and E = respectively; they are 
depicted in Fig. 7b, and 7c. Note that their kinetic en- 
ergy is Ek =0; their only energy comes from the chemical 
potential. 



B. Second order perturbation theory and 
renormalized hopping and tilt 

Having diagonalized the bond Hi we now need to in- 
corporate the rest of the Hamiltonian as a perturbation. 
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W- W- Wj 



FIG. 7: Eigenstates of the bond i. (a) The lowest eigenstate of 
the Hamiltonian TLi from Eq. (31) is |— ) (Eq. 35). Sites i and 
i + 1 share a single vortex. In the imaginary time description, 
it is as though the vortex tunnels back and forth between the 
two columnar pins. The energy of this state is — /j, — Wi. (b) 
An excited states in which both pins are empty, with energy 
eigenvalue zero, (c) An excited eigenstate of Tit with energy 
—2p,\ both pins are occupied by vortices. 




FIG. 8: Real-space RG decimation step, where we allow for 
different tilt fields hj on neighboring bonds. In each step 
of the renormalization we find the strongest bond, Wi — Q. 
We diagonalize it and minimize its energy by allowing a single 
vortex to fluctuate between the two columnar pins, i and i+1. 
A vortex in site i — 1 may then fluctuate into site i, raise the 
energy of sites i and i + 1, and make the vortex in site i + 1 
relax by fluctuating into site i+2. Thus we can eliminate sites 
i and i + by including the effective hopping between sites 
i — 1 and i + 2. The effective transverse field (tilt) between 
i — 1 and i + 2 becomes the sum of the tilts connecting these 
two sites. 



To second order in the hopping matrix elements, we need 
to concentrate only on the two neighboring bonds: 



V 



m-i (bj^bie h + bTbi_ie 



-w. l+1 [h\ +l h l+2 e 



-'j-ht h i if 



(37) 



The second order Hamiltonian is then: 

y(2) = (-|V|0,0)(0,0|V|-) + (-|V|l,l)(l,l|V|-) 



(38) 



h\_ x h i+2 e- 5h + bt^bi-re 



■Mi 



where, as in Eq. (36), |0,0) is a state with both sites i 
and i + 1 empty, and |1, 1) is a state with both sites full. 
Quantum fluctuations thus give rise to a new coupling 
between sites i—1 and i + 2. This renormalization process 
is depicted in Fig. 8, and is discussed below. 

Eq. (38) is quite remarkable. It tells us that after 
localizing a vortex between sites i and we can forget 
about these two sites. Except for a constant contribution, 
the only change we need to make to the Hamiltonian is to 
add a new hopping term that allows vortices to hop over 
the occupied pair, from i — 1 to i + 2. In addition, we see 
that this hopping is associated with the composite tilt of 
the three bonds linking i — 1 and i + 2. In general, upon 
allowing for different tilt fields on neighboring bonds, the 
effective strength and effective tilt of the renormalized 
bond are respectively 



eff 



-2 — n 



h e ff — h 

n i-l,i+2 - n i-l 



(39) 



where we introduced the notation = max{wi} (also 
see Fig. 8). Eqs. (39) are the RG rules for the RSRG as 
applied to the vortex hopping problem. 

As it turns out, Eqs. (39) coincide with the RG flow 
equations for the fermion hopping problem (or XX spin 
chain problem) in the Hermitian case h = 0. The only 
difference is that the tilt hi replaces the length of bonds, 
as one might expect in the case of localized states. 2 In or- 
der to solve the flow equations we can thus use techniques 
known from the Hermitian case. 

Two important points should be made here: (1) The 
flows of the hi's and w^s are independent. (2) No 
new terms are produced. In particular, we do not pro- 
duce randomness in the chemical potential. In the pure 
model it is the particle-hole symmetry that prevented us 
from having new chemical-potential terms. But Eq. (7) 
with nonzero tilt doesn't have a particle-hole symmetry, 
or more precisely, it is symmetric under a particle-hole 
transformation accompanied by a changing the tilt direc- 
tion: 

b -► h\ bt -> b h-> -h (40) 



Terms like 



H+l 



b f be' 1 + bbV* 



are allowed in the presence of the modified p-h symmetry 
(Eq. 40), and they would create a bias towards vortices 
(or holes). Fortunately, as discussed above, they are not 
produced. 
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C. Flow equations for the distribution of hi and w» 

In the previous section we derived the RG rules that 
govern the effective tilt and hopping matrix elements. By 
using these RG rules we can gradually decrease the en- 
ergy scale of the Hamiltonian while determining the high 
energy parts of the ground state. When the energy scale 
reaches zero, the Hamiltonian that we started with is ef- 
fectively diagonalizcd. The RG rules on their own, how- 
ever, do not tell us much about the eigenvalue spectrum. 
To obtain useful information, we need to convert the RG 
rules into differential flow equations for the distributions 
of bond strength and tilt. As pointed out before, this 
has been essentially done already in the corresponding 
Hermitian problem. 

Eq. (39) provides RG rules that are identical to the 
case of the random XX Hamiltonian. 13 In this case we 
rcnormalize the bond strength and the length of bonds. 
Eqs. (39) suggest that for random vortex hopping in the 
non-Hermitian problem, the length of a bond is replaced 
by the bond's tilt. This observation allows us to use es- 
sentially all results derived for the XX model in Ref. 13 
and apply them in this problem. In this section we re- 
view these results as they apply to the vortex hopping 
problem. 

We need the flow equation for the joint probability 



distribution p(w, h) of the hopping strength and tilt on 
a bond with respect to the gradually decreasing energy 
scale fi. As in the theory of random spin chains, we 
introduce a logarithmic variable Q to replace the bond's 
strength, and a renormalization-group flow parameter V: 



»' Wi ' 

r — In ill! 
1 111 SI ' 



(41) 



where fio is the initial energy scale (largest matrix ele- 
ment) in the problem. The logarithmic energy scale T 
is the RG flow parameter, and it increases as the energy 
scale is decreased. Note also that in every step of the 
RG we renormalize away bonds with ( — > 0. In the con- 
tinuum limit, one can verify that to generate the flow 
r — > T + dT we renormalize all bonds with < ( < dT. 
The logic behind the definition of ( in Eq. (41) becomes 
clear upon considering the first RG rule in terms of the 
new variables. Eq. (39) takes the form 



C eff - C- 



(42) 



Upon transforming p(w, h) into the new variables, we 
obtain the joint probability distribution P((, h). The 
flow equation for P((, h) is given by 13 



J 



dP(c,h) _ dP(c,h) 



dT 



+ 



J J dCid^dhodh^SiC - Ci - &)5(h -ho -hi- h 2 )P(0, h )P(( u h 1 )P(C2,h 2 ). (43) 



Eq, (43) greatly simplifies if we take the Laplace trans- 
form with respect to h. If y is the Laplace transform 
variable conjugate to h, we obtain 

dP((.,y) _ dP(C,y) 

+ JfdC 1 d( 2 6(( - Ci - C 2 )P(0, y)P(C 1 ,y)P(C2,y)- 

(44) 

The derivation of this equation is given in detail in Ref. 
13. 



D. Universal fixed point distributions of w and h, 
length energy scaling, and density of states 

Remarkably, Eq. (44) has a scaling solution which is an 
attractor to essentially all boundary conditions. 14 This 
solution is 

p(C, y) = t^-tt cxp (-Vv coth (v» (r + r ) 0) (45) 

1 + 1 

where Tq is a non-universal constant of integration. This 
solution, once inverse-Laplace transformed, gives the 



joint probability distribution of bond strengths and ef- 
fective tilt. Given some initial Hamiltonian and tilt, at 
sufficiently low energy scales (i.e. in the limit of large 
system size for the vortex problem) the distribution of 
effective couplings will be given by Eq. (45). 

From this solution we find that the average logarithmic 
coupling obeys: 

C = r + L , (46) 

where over-bar denotes disorder averaging. Also, Eq. 
(46) is obtained from the marginal distribution for the 
logarithmic couplings, i.e., Eq. (45) with y — > 0. In ad- 
dition we infer that the Laplace parameter conjugate to 
the tilt scales as 

y~ (r + r )- 2 , 

and therefore the tilt field scales as 

7^77- (r + r ) 2 . (47) 

From the above scaling we can infer the average length of 
a bond after the RG reached the logarithmic energy scale 
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T. This scaling is identical to the tilt-energy scaling, as 
mentioned above: 



/ = / (r + r ) 2 



(48) 



The concomitant relation to the length-energy scaling, is 
the relation between the fraction of unoccupied pinning 
sites, /, and the energy. This relation has to be 



/ = 



fo 



(r + r ) 



2 ■ 



(49) 



where r and f are non-universal constants that depend 
on the particular realization of the disorder. 

A formula used by Shnerb and Nelson, 8 ' 9 and Brouwer 
et al. 6 will provide us with a connection between the spec- 
trum of the Hermitian Hamiltonian (with zero tilt) and 
the spectrum of the tilted Hamiltonian, see Sec. VI B. In 
order to use this formula to find information about de- 
localized vortex states, however, we will need knowledge 
of the density of states in the zero-tilt limit. 

The density of states (DOS) can also be derived from 
Eq. (49). The number of states per length in the energy 
interval Q -> - dfl, or T -> T + dT (with ft = fi 7 e- r ) 
is just the number of bonds that get decimated at that 
energy scale. This, in turn, is given by: 

/■P(C = 0,y = 0)dT. 

Converting this into energy terms, using Eqs. (45) and 
(49), we obtain 



9(e) = 



fo 



1 



\n 4 (fij/e) e 



(50) 



where the energy e is measured from the center of the 
band, and 



0/ = O e °. 



(51) 



We expect Eq. (50) to be asymptotically exact for small 
e. 

In order to check Eqs. (46) and (49), we carried out 
the real-space RG procedure numerically on chains with 
5 x 10 6 sites for two types of distributions - block distri- 
bution and displaced power law distribution. In Figs. 9 
and 10 we plot the average logarithmic coupling and the 
density of free sites vs. the logarithmic flow parameter, 
r, for the two distributions types (further described in 
the figure captions). From these plots we also infer the 
values of the non- universal constants fo and r , which 
will be used later to compare results for the vortex local- 
ization problem. 



E. Applicability of the RSRG with nonzero tilt 



no longer write a vortex wave function which is local- 
ized between two sites - this becomes impossible as soon 
as the chemical potential reaches the mobility edge, and 
the states become extended. We can nevertheless take 
advantage of the real-space RG as long as we stop deci- 
mating when the energy scale reaches the mobility edge. 
In addition, the RSRG gives a reliable expression for the 
density-of-states of the Hamiltonian Eq. (7) with zero- 
tilt. 



V. CONNECTION WITH VORTEX PHYSICS 

In the previous section we derived the RSRG for the 
vortex-hopping imaginary time Hamiltonian. In this sec- 
tion we will clarify the relation of the results obtained 
above to the physics of vortex lines with columnar pins 
in planar superconductors. 

The vortex hopping Hamiltonian yields an energy spec- 
trum, with filling factor determined by a chemical poten- 
tial, \x. The RSRG outlined above is used to determine 
the wave function of vortices in a given chemical poten- 
tial; in each step of the RG, the strongest bond is diago- 
nalized, and if the bond-energy of a vortex placed on the 
bond is lower than the chemical potential, 



E- = -w - ii < 0, 



(52) 



then a vortex is localized on that bond. If the strongest 
bond in the renormalized chain is too weak to obey (52), 
then the RG is stopped and we have found the ground 
state of the zero-tilt model. 

For flux lines, however, we often want to know the 
longitudinal magnetic field, proportional to the vortex 
density n. We need to find out how n and fj, are related. 
This is done using Eq. (49). The longitudinal vortex 
density vs. the fraction of unoccupied sites / is 



f) = \ 



fo 



(53) 



i.e. the longitudinal vortex density is half the density 
of filled pinning sites, since each vortex is localized on 
bonds, and thus on two sites. The number of unoccu- 
pied sites / as a function of energy scale is taken from 
Eq. (49). Now, the condition (52) gets translated to 
(assuming negative /x; for positive fx we can carry out a 
particle-hole transformation) 



T = In- 



ImI 



Upon inverting Eq. (53), we thus obtain 



= — f2oe r ° exp 




(54) 



(55) 



The application of the RSRG breaks down as soon as 
we reach energies at the mobility edge, because we can 



Note that flo is the initial energy scale, and fo and To are 
constants that are determined by the initial distribution 
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FIG. 9: Numerical real-space RG for an initial block distribution, i.e., initial hoppings uniformly distributed between d and 
1 + d. Each graph represents averaging over 20 realizations, (a) The average logarithmic coupling, ( vs. the logarithmic RG 
flow parameter, F. The fitted slope in the two curves is 0.99. From the intercept we deduce Fo for the two curves: To = 0.31 
(d — 0.1) and To = 0.19 (d = 0.3). These plots should be compared with the linear dependencies predicted by Eq. (46). (b) 
The inverse density of free sites 1// vs. (r + To) 2 , with To obtained from the graphs in (a). As expected from Eq. (49), the 
intercepts of the curves with the y axis in the two plots are negligible compared to 1//: 2.6 (d = 0.1) and —2.4 (d — 0.3). From 
the slope of the two curves we find: fo = 0.37 (d = 0.1) and fo = 0.16 (d = 0.3). 



a=0.5 






FIG. 10: Numerical real-space RG for an initial 'displaced power- law' distribution, with initial hoppings given by w = 0.01 + w r 
and < w r < 1 distributed as P[to r ] = (1 — a)^ja- The distribution of w r corresponds to a hopping strength which is 

exponentially suppressed with the distance between pins: w r oc e~ x with x distributed as P[x] oc e - ' 1- ™' 1 . Each graph 
represents averaging over 20 realizations, (a) The average logarithmic coupling £ vs. the logarithmic RG flow parameter F. 
The fitted slope in the two curves is 0.99. From the intercept we deduce Fo for the two curves: To = 0.71 (a = 0.5) and 
To = —0.08 (a = 0.8). As before, these plots should be compared with Eq. (46). (b) The inverse density of free sites 1// vs. 
(r + To) 2 , with Fo obtained from the graphs in a. As expected from Eq. (49), the intercepts of the curves with the y axis in 
the two plots are negligible compared to 1//: —0.29 (a — 0.5) and —0.6 (a = 0.8). From the slope of the two curves we find: 
fo = 1.85 (a = 0.5) and fo = 2.6 (a = 0.8). 
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n 

FIG. 11: Semi-log plots of the vortex chemical potential, /i, 
vs. the longitudinal vortex density, n, for a box disorder dis- 
tribution and two values of d. The solid lines are Eq. (55) 
with To and fo given in the caption of Fig. 9. Insets show 
the linear-linear plots. 



of hoppings, and which can be calculated as in Fig. 9. 
Eq. (55) is tested in Figs. 11 and 12 for two forms of 
disorder. 

In the following, we will find the mobility edge, /i m , of 
the vortex- hopping model with tilt, using the RSRG. We 
will then use the above relations to convert the results to 
the physical variable n. 

All results will depend on the initial conditions through 
the non- universal quantities £1/ and / . Nevertheless 
these results are valuable since our results will be able 
to describe the entire phenomenology of the depinning 
transition in terms of these two constants, which encode 
the relevant part of the initial distribution. These con- 
stants can be evaluated for a given distribution either by 
carrying out the RSRG numerically, or by fitting curves 
to partial results obtained in other methods, for instance, 
by comparing the density of states, Eq .(50), to a spec- 
trum obtained numerically. 

In general, when comparing our results to spectra ob- 
tained numerically, it is easier to use formulas with ji in 
them rather than the physical n. 



n 

FIG. 12: Semi-log plot of the vortex chemical potential, fj,, vs. 
the longitudinal vortex density, n, for two displaced power- 
law distributions as in Fig. 10. The solid lines are Eq. (55) 
with To and fo given in the caption of Fig. 10. Insets show 
the linear-linear plots. 



VI. APPLICATION OF THE RSRG RESULTS 
TO NON-HERMITIAN DELOCALIZATION 

In the above we derived many characteristics of the 
vortex hopping problem while assuming that all vortices 
are localized. The most interesting aspects of the prob- 
lem, however, arise in relation to the derealization tran- 
sition of vortices with sufficient tilt. 

The physical quantities of interest are essentially all re- 
lated to the energy-spectrum of the problem. Most par- 
ticularly, given a tilt - an angle relative to the columnar 
defects with which we apply a magnetic field - how much 
magnetic field can we apply such that all vortex lines re- 
main localized in configurations parallel to the columnar 
pins? This question is answered by finding the mobility 
edge of the non-Hermitian boson hopping model. Once 
we know the effective chemical potential in which the 
vortex bosons become depinned, the density of bosons 
at this chemical potential determines the critical field for 
destruction of the transverse Meissner effect. 

Another interesting quantity is the transverse magnetic 
field as a function of tilt above the mobility threshold for 
a given longitudinal field. The transverse field is propor- 
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tional to the boson current [see Eq. (8)]. For tilts slightly 
above the mobility threshold h c (i.e., a chemical poten- 
tial slightly above the mobility edge), the boson current 
is simply given as (Eq. A4): 

J = — (/j — fj, c ) tan 9 

7T 

where 9 is the angle of ascent of the spectrum in the 
complex energy plane (see Fig. 16). The angle 9 is a 
universal function of the external magnetic field, or the 
chemical potential. 

In this section we will analytically derive these quanti- 
ties for the random hopping model. First, we will find the 
mobility edge as a function of tilt, or, equivalently, the 
tilt threshold as a function of external field. In order to 
obtain the transverse magnetization, we will need to use 
a formula appearing in Refs. 6 and 8,9 , which connects 
the spectrum of the Hermitian boson hopping Hamilto- 
nian with the spectrum of the non-Hcrmitian problem 
and the tilt. We will introduce this formula, and use it 
to derive the quantities of interest. 



Mobility edge of the random hopping 
Hamiltonian 



By using the RSRG, we demonstrated in Sec. IV that 
when the tilt is zero, vortices are localized between pairs 
of pinning sites. We also showed that as we reduce the 
energy scale of the Hamiltonian, and place vortices onto 
pairs of pinning sites that are increasingly far apart, the 
effective tilt, h e s per site also increases. For any given 
initial tilt, at some energy scale there will no longer be 
a localized vortex solution. This energy is the mobility 
edge, n(h). Alternatively, there will be a critical longitu- 
dinal vortex density, which we denote n(h), above which 
vortices are no longer restricted to the pinning sites. In 
this section we find the relation between the tilt and the 
critical longitudinal vortex density. 



U 




Tilt direction 



FIG. 13: Simplified model with a single strong pin. We use 
this model to determine the mobility edge for a given tilt. 
The model consists of a uniform pinning lattice with hopping 
w, periodic boundary conditions, and one pair of strongly in- 
teracting sites with hopping wo. The strong bond localizes 
a vortex in it as long as the tilt is sufficiently weak. The 
wave function is schematically plotted, showing that the its 
decay is asymmetric. Strong enough tilt will prevent the wave 
function from decaying on its right side, signaling a dereal- 
ization transition. This problem with non-hard-core vortices 
was considered in Ref. 7. 



The eigenfunctions of Hamiltonian (56) can be solved 
exactly. In particular, the localized wave function at the 
bottom of the band reads, up to normalization, 

e M"-V2) b t + £ e—<-V»>bt) |0) . 



V n<0 



n>0 



(57) 

The left and right decay coefficients, kl and kr, are given 
by: 



K R a = In 2& - h, 
n L a = In ^ + /i, 



(58) 



where a is the lattice constant. 

A localized solution of the form Eq. (57) can only exist 
as long as both kr and kl are positive; otherwise we 
cannot accommodate the periodic boundary conditions 
(as the system is a ring). Hence, the threshold tilt is 
given by the condition that kr = 0, or 



1. Single strong bond 



In — . 

w 



(59) 



We start our investigation by solving a simplified prob- 
lem that will give us the correct solution. We will then 
proceed to justify it for the random hopping Hamilto- 
nian. Consider a ring of columnar defects, where there 
is one strong bond, Wo, and all other bonds are w <C Wo 
(Fig. 13): 



H = -E» (hlb l+1 e- h + bl +1 b ie h ) 



i#0 
-Wq 



(bJbie-^ + blboe' 



(56) 



The strong bond can localize a single vortex; we will 
determine at what tilt h this state delocalizes. 



2. Generalization for the random-hopping model 

Let us now generalize this result to the case of ran- 
dom hopping. At every stage of the RG we assume 
that the strongest bond localizes a vortex, with energy 
E = Q = £! e~ r . If we assume that all other bonds 
are the same, the problem reduces to the localizing bond 
problem above. However, we need a generalized condi- 
tion that takes into account the randomness of the weak 
bonds. In this case the natural generalization of condi- 
tion (59) would be: 



t off 



w w 



To. 



(60) 
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where the last equality is due to Eq. (46). The reason 
for writing h en ? is to remind ourselves that the tilt h we 
consider in this section is the renormalized h. In the 
next section we connect the results derived above with 
the bare parameters of the vortex-pinning problem. 

So far Eq. (60) is just a physically motivated guess. 
There are, however, several more rigorous ways to derive 
it. The most straightforward proof is obtained via per- 
turbation theory. Consider the Hamiltonian Eq. (31), 
and assume that coupling w a is by far the strongest in 



the chain. From Eq. (34) we know what the zeroth order 
solution for the wave function is: 

|_) = -L(e-^bS + e fe/2 b t 1 )|0) (61) 

From perturbation theory, we can easily see that the wave 
function at site |m), to lowest non- vanishing order in the 
Wi's reads (assume m > 0) 



(m|-) 
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\i=2 



(iKe fc btbi_i|i-l)\ (2|tuie h b5bi|-) _ e 



oh/2 



W 



W 



V2 



n 

i=l 



e'-Wi 



w 



(62) 



By requiring that the product at the end of Eq. (62) 
remains finite as m — > oo, we obtain condition (60). An- 
other derivation is supplied in Appendix B. 



3. Critical tilt and field 

As explained in Sec. (IV) above, the RSRG eliminates 
the strongest hoppings in the Hamiltonian (7) and grad- 
ually populates localized vortex states, until the energy 
scale of the renormalized Hamiltonian reaches the chem- 
ical potential: 

Q = max{w,} = — fi, 

i 

where the tilde indicates that the maximum is evalu- 
ated over the set of remaining (and renormalized) bond 
strengths. The chemical potential, /i, in this case is a 
tuning parameter that controls the vortex filling factor 
of the lattice, and hence the longitudinal vortex density. 
In terms of the RG variables, £1 = Q e~ r , and therefore 
we can define 



hA 

-n 



(63) 



From the RG procedure, and in particular from Eq. 
(47), we know that 



h e a ~ h (r M + r ) /, to- 



rn 



where To and fo were defined in Eqs. (46) and (49) re- 
spectively, and ho is the unrenormalizcd tilt per lattice 
site. 

In order to find a relationship between the mobil- 
ity edge of the vortex hopping problem and the un- 
renormalized tilt h, we substitute the relation h°fj — 

h c (r M + T ) 2 /fo into Eq. (60) and obtain 



h c 



fo 



(65) 



The physical quantity in the vortex pinning problem 
that the parameters /x and T M determine is the vortex 
density, n. As a function of T M , we showed in Sec. (V) 
that 



n = - 



fo 



2(r„ + r r 



and therefore the connection between the critical tilt and 
the longitudinal field is 



h c = Vfo (1 - 2n) 



(66) 



This is the main result of this section. At a given field 
n, as we increase the bare tilt, when the tilt reaches h c , 
vortices get delocalized and transverse flux ensues. The 
critical tilt vanishes asm 0.5, i.e., as vortex states fill 
the lattice upto the middle of the band. This is since the 
localization length of the tilt-free vortex states diverges as 
the chemical potential approaches zero. We demonstrate 
the validity of Eq. (66) using the numerical real-space RG 
in Figs. 14 and 15. Note that the RSRG approach is most 
accurate near the middle of the band, where n <~ 0.5, in 
the region where the density of states diverges as in Eq. 
(50), also see Sec. VII. 



4- Vortex localization length 

Another quantity of interest is the vortex localization 
length. According to Rcf. 2, this length, is de- 
fined as the decay distance of the ket describing the least 
localized vortex in the tilt direction. From Eq. (58) 
and its generalization according to Eq. (60) we see that 
£_l = — ~ n 1 . , where the tilde indicates that the 



«H In kl-h eff ' 

length is a renormalized length, and not bare length. It 
follows that 



1 



Tn he 



(67) 
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FIG. 14: Critical tilt, h c , vs. longitudinal vortex density, n, 
for the two box-distributions from Fig. 9a. The dots are the 
result of the numerical implementation of the RG, where the 
critical tilt is the sum of £ divided by the full (bare) length 
of the chain. The solid line is a plot of Eq. (66). 




n 



FIG. 15: Critical tilt, h c , vs. longitudinal vortex density, n, 
for the two displaced power-law distributions from Fig. 10. 
The dots are the result of the numerical implementation of 
the RG, where the critical tilt is the sum of £ divided by the 
full (bare) length of the chain. The solid line is a plot of Eq. 
(66). 



To convert the renormalized length to physical length, 
we need to divide by the fraction of empty pinning sites, 
given by Eq. (49). In terms of physical (bare) length 
scales, the localization length is 



B. Non-Hermitian spectral formula 

A formula in Ref. 8,9 (see also Ref. 6) relates the 
imaginary part of the eigenenergies of delocalized states 
with the spectrum of the Hermitian hopping problem, in 
the absence of a tilt. This connection, along with our de- 
tailed knowledge of the random-singlet phase allows us to 
probe even the delocalized states, even though the RSRG 
does not strictly apply for these states. This remarkable 
state of affairs stems from the analytic properties for a 
next-neighbor hopping problem. 

The non-Hermitian spectral formula is discussed exten- 
sively in Ref. 8,9, but we quote it here for completeness: 

JV N 

H(E-e i ) = 2 [cosh(hN) - 1] JJ (- Wi ) , (69) 

i=l i=l 



where E is an eigenvalue of the non-Hermitian vortex 
hopping problem with tilt h per site. There are a total of 
N sites arranged in a ring. The e^'s are the eigenvalues 
of the zero-tilt Hermitian problem. This is a complex 
equation, and we will use it primarily to hnd a relation 
between the imaginary part of E and the the excess tilt 
h — h c . In principle, we can use the above formula to 
find all eigenvalues of the non-Hermitian problem, since 
we know the spectrum of the zero-tilt problem from the 
RSRG. This strategy can be implemented numerically. 

A particularly useful form of the non-Hermitian spec- 
tral formula is as follows: Let us assume that N and Nh 
are large, and take the logarithm of Eq. (69). We obtain 

JV JV 

J2^(E-ei) = N\h\ + J2^(wi)+iTTN. (70) 

i=l i=l 

In Appendix B we show how the RSRG result Eq. (60) 
can be deduced from the discrete non-Hermitian spectral 
formula. There it is also shown that we can, in fact, 
apply the RSRG directly to formula (70). The advantage 
of doing so is that in advanced stages of the RSRG the 
distribution of Wi is known from the fixed point solution, 
Eq. (45). 
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FIG. 16: (a) Typical spectrum of the random vortex-hopping 
Hamiltonian with non-zero tilt. The energy eigenvalues lie 
on the bold lines. Localized states lie on the real axis, while 
delocalized states have an imaginary component. Qo is the 
strongest bond in the bare Hamiltonian, and is of the order of 
Qi, which is the parameter used to define T - the logarithmic 
energy scale. 9 is the angle of approach of the complex branch 
of the spectrum, and is a critical property of the depinning 
transition, (b) To calculate 9 we concentrate on the wedge at 
the mobility edge. The bold line is the spectrum at tilt h + dh 
and the dashed is the complex spectrum at tilt h. 



An important formulation of Eq. (70) is obtained by- 
taking its real part: 



N 

E 



i In ({RcE - e t ) 2 + {ImE) 2 ) = N\h\ 



N 

E 



In (uii) 



(71) 

Eq. (71) specifies a curve in the complex E-planc on 
which the eigenvalues of the vortex hopping problem lie. 
In the absence of the imaginary part of this equation 
it is impossible to determine the location of individual 
eigenvalues. But if we are interested in the shape of the 
curve, Eq. (71) is all that is necessary. 

In the following we will make use of the above formulas 
in the context of the random hopping problem to obtain 
analytical expressions that describe the depinning tran- 
sition and the spectrum of the delocalized states. 



to the current of delocalized 'bosons' (i.e., the number of 
tilted vortices), as is discussed in the next section. 

The calculation of the angle of approach is somewhat 
subtle, and we refer to Fig. 16b to aid the discussion. 
The figure shows the line on which the spectrum lies when 
the tilt is h + dh in bold. The dashed bold lines indicate 
the line of spectrum for the lower tilt h. The angle of 
approach is given by the ratio of the segments: t&nO = 
AB/AC. 

We first determine the vertical distance between points 
A and B in the figure. For this purpose we differenti- 
ate the real part of the non-Hermitian spectral formula, 
Eq. (71) with respect to ImE while keeping RcE con- 
stant. To avoid confusion we use a representation in 
terms of differences, taking the difference between the 
expression with ImE = c and ImE = c + de. Upon 
passing to a continuum representation and introducing a 
density of states g(e), we obtain 



de-g{e) 



c ■ de 



(RcE 



e 4 ) 2 +c 2 



= dh. 



(72) 



-n 

In the limit of dh — ► 0, the integrand on the left is nothing 
but g(e) ■ de ■ 7r<5( ReB _ e ), and we obtain (where we drop 
the Re sign) irg(E)de = dh and hence find that AB in 
Fig. 16b is 

,3 fij 



de = dh ■ 



\E\ In 6 



(73) 



where g(E) is taken from Eq. (50). In Eq. (73) E desig- 
nates the energy of the mobility edge, so we substitute 

E = Llm(h). 

The second part of the calculation focuses on AC. 
This is the distance the mobility edge travels as the tilt 
changes from h to h + dh. But from Eq. (65) we can find 
this quantity directly: 

dh _ f a fo 



I In 2 FT 
1 M 



d^ m | M | (T„ + To) I Mr, 
The angle of approach is now readily found: 



tan 6 = 



de 



dh 



de 
dh 



7T 



ft/ 



(74) 



(75) 



dn m d^„ 

Again it is a universal form, which depends on the initial 
distribution only through a single parameter, ftj. Note 
that as the tilt h approaches zero, and the mobility edge 
approaches the 'middle of the band' (i.e., [i m — ► 0), the 
angle 9 approaches n/2. 



D. Vortex current near the mobility edge 

C. Angle of approach at the mobility edge 

The imaginary vortex current near the onset is in gen- 
An interesting quantity, which enters the critical prop- eral given by Eq. (A4) in App. A: 
crties of the depinning transition, is the angle of approach ^ ^ ^ 

of the delocalized branch of the spectrum relative to the Jtotai — 2— Re — = — (ImE^ — ImE^ ) , (76) 

real axis, as shown in Fig 16. This quantity is related ^ kn ®^ tt 
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where L is the number of sites in the lattice. In partic- 
ular, we can apply this formula to the random hopping 
problem. In Eq. (75) the angle of approach near the 
delocalization transition was found (also see Fig. 16): 



100 E 



dhnE 
de 

where we define e as: 



tan# 



1 . Hi 
-In- : 

7T | Mm | 



e = ReE, 

and Hm is the chemical potential at the mobility edge. 
Now, near the onset of transverse flux penetration we 
can write: 



J total 

i tan#(/x 



) = 



= -4 In 



1^1 (M- Mm). 



(77) 



This result can also be expressed in terms of the dif- 
ference between the longitudinal vortex density and its 
critical value. The number of extra vortices per length is 
just n — n c . Therefore the transverse vortex density, i.e., 
the imaginary 'boson' current is also: 

1 T de tnnO ,_ . 

Jtotal ~ T Im ^7 \ n _ n c) ~ i n - n cj 7 T ('8) 

L Ok tt ■ g(fi m ) 

Close to n c , we may assume that the real parts of the en- 
ergy of eigenstates near the mobility edge do not change 
as they undergo the delocalization transition. Therefore 
we can use Eq. (50) for the DOS. Upon combining Eq. 
(50) with Eqs. (55), (65), and (66) we can write result 
(78) in terms of the longitudinal vortex density only: 



Jt 



total 



(n 



i 



(1 - 2nf 



= -\//o/(l-2n) 



or in terms of the tilt only: 



Jtotal ~ (h — h c ) 



7T 2 h 3 ' 



-fo/h 



(79) 



(80) 



Eqs. (79) and (80) are the main results of this section. 

One can actually use the above reasoning to obtain 
not only the total current, but the current of individual 
states near the mobility edge. Consider the state E n 
which can be associated with a wave vector k n - As we 
change k = k n to k n + =jr-, we move from the eigenvalue 
E n to E n+ \. Thus 



Re(dE) 



tn+i 



L-g(e) 



(81) 



where we make use of Eq. (50), and again set e = Re(E). 
The factor 2 in the numerator of Eq. (A3) is due to the 
real spectrum splitting into two branches, and we choose 
a convention in which both branches of the spectrum 
contribute to the DOS. Now we can also substitute dk = 
=jr- t and we obtain for the current of this state 



j - iRe^ 



L-g(e) 2tt 



dh 

tan 9 



(82) 
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FIG. 17: Density of states of the random hopping problem 
with zero tilt for 2000 lattice sites, averaged over 50 dis- 
order realizations. The hopping parameters Wi are equally 
distributed in [0.2,0.7]. At low energies, good agreement is 
obtained with the analytic RSRG result Eq. (50). 



where we used tan 9 



dImE 

de 



The result (82) can be 
written in terms of the longitudinal vortex density, n, 
instead of in terms of e: 



1 



Ltt 2 (1 - 2n) 



(83) 



Eq. (82) can also be expressed in terms of the tilt h, 
using Eq. (66): 



j ~ fo^ 1 e -Jb/fc 
11 ~ Ltt 2 h 4 



(84) 



VII. 



EXACT DIAGONALIZATION OF FINITE 
SYSTEMS 



In this section we analyze the free-fermion Hamiltonian 
(13) by exact numerical diagonalization of finite-size sys- 
tems with up to N = 4000 lattice sites. While this study 
is mostly focused on the case of random hopping, we 
also present some results for on-site disorder. We employ 
periodic boundary conditions in order to model qualita- 
tively a realistic superconducting slab where vortices can 
enter and exit at the edges. Note that for open bound- 
ary conditions (i.e., zero hopping matrix elements at the 
edges), the magnetic field tilt could be gauged away and 
would not induce any transverse vortex density. Since all 
quantities are random, we average over a finite number of 
disorder realizations, Nr — 20 . . . 100. Unless indicated 
otherwise, the random hopping matrix elements Wi are 
drawn from a flat distribution in [c, c + 0.5], where the 
lower cutoff c = 0.2 is used to prevent singular couplings 
close to zero, which would lead to almost decoupled sub- 
systems. For diagonal disorder, we use on-site energies 
ti which are evenly distributed in [—0.5,0.5]. 



19 



100 



10 



- 1 1 1 






o L=500 






1000 






2000 






-- (1/0.65) (ln(0.8/lel)) 2 







IK 



le-05 0.0001 0.001 0.01 0.1 



1000 



100 



10 




L=500 

1000 

2000 



0.1)01 



0.01 



0.1 



FIG. 18: Top: localization length of the random hopping 
problem with zero tilt for different system sizes L, aver- 
aged over 50 disorder realizations. As in Fig. 17, the hop- 
ping matrix elements Wi are equally distributed in [0.2,0.7]. 
The low-energy behavior agrees very well with the result 
£ ~ ± ln 2 (f2//|e|) predicted by the RSRG. Bottom: local- 
ization length for onsite disorder with site energies et equally 
distributed in [—0.5,0.5], averaged over 10 realizations. £(e) 
clearly stays finite as e — > 0. 



We first study the case of vanishing tilt h = 0, in order 
to illustrate the qualitatively different physics of bond 
and on-site disorder. For bond disorder, we have calcu- 
lated the density of states from the real part of the single- 
particle spectrum, averaging over Nr — 50 realizations 
of the disorder. The results are shown in Fig. 17, where 
it is evident that the DOS diverges at the center of the 
band, as predicted by the RSRG. We obtain a good fit 
to the theoretical result in Eq. (50). 

A similar behavior is obtained for the localization 
length £, which we define as 



L 

E 



*i\ 2 \*i\*dc(i,3) 



(85) 




0.0001 



FIG. 19: Spectrum of single-particle states for different values 
of the tilt, with hopping disorder as in Fig. 17, L — 2000 
lattice sites and averaging over 100 disorder realizations, on 
a linear scale (top) and in a lin-log plot (bottom). 
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FIG. 20: Mobility edge vs. tilt for different values of the tilt, 
with hopping disorder as in Fig. 17, L — 2000 lattice sites and 
averaging over 100 disorder realizations. 



20 




FIG. 21: Imaginary current for random hopping distributed as in Fig. 17, as a function of tilt h and vortex density n. Results 
are given for L — 1000 sites, with a disorder average performed over 20 realizations. 



0-02 




FIG. 22: Onset of imaginary current (i.e., transverse magnetic 
flux) as a function of tilt, for random hopping as in Fig. 17 
and vortex density n = 0.25. Data are averaged over 20 real- 
izations of disorder. Apart from finite-size rounding the onset 
is linear, as predicted by Eq. (80). 

where 

dci iJ) = L skl (^A^ ( 86) 

is the chord distance on the ring and \&j is the wave 
function at site i of the single-particle state under con- 
sideration. Results are shown in Fig. 18. For random 
hopping we find that the localization length diverges in 



the middle of the band according to 



as predicted by the RSRG [Eq. (48)]. In particular, 
there is always a vortex state at zero energy which is 
delocalized across the whole system. On the other hand, 
for on-site disorder £ clearly stays finite, as shown in the 
right graph of Fig. 18. As a result, for bond disorder an 
arbitrarily small but finite tilt h leads to extended states 
at the center of the band, in contrast to the case of on- 
site disorder where this only happens above a critical tilt 
h > h c . 2 [Note that the localization length discussed here 
is the distance between the two pinning sites that share 
a vortex, as opposed to £j_ of Sec. VI A, which describes 
the wondering of the vortex away from the pinning sites 
due to the tilt.] 

We now present results for finite tilt h > 0. In this 
case the Hamiltonian is non-Hermitian, with complex 
eigenvalues and left/right eigenstates which are no longer 
equal. In Fig. 19 we show numerical results for the 
spectrum of the random hopping problem. For small 
tilt almost all states are localized and the corresponding 
eigenenergies real. As the tilt increases, the two mobil- 
ity edges ±/i c move towards the band edges and define a 
growing region of extended states. In Fig. 20, the depen- 
dence of fj, c on the tilt is shown. The rapid vanishing of 
fi c at small h is consistent with Eq. (55) using Eq. (66). 
(two fits are shown), although a quantitative comparison 
is difficult due to the finite system size. 
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At finite tilt, the most interesting physical quantity for 
the vortex problem is the transverse flux, corresponding 
to the imaginary current J. Following Eq. (14), we obtain 
the current as the derivative of the ground-state energy 
Eg with respect to tilt: 



T ( ^ ldE 

J={ - l) L^h 



(88) 



The resulting current for hopping disorder, as a func- 
tion of particle (vortex) density n and tilt h, is shown 
in Fig. 21. Clearly, for weak tilt and few vortices, the 
current vanishes since the mobility edge is close to the 
center of the band and all the occupied states are local- 
ized. This is a manifestation of the transverse Meissncr 
effect: In the presence of disorder due to columnar pins, a 
weak transverse magnetic field does not induce a trans- 
verse magnetic flux, because the flux lines are pinned 
by the defects. 1 With increasing filling and tilt, a well- 
defined transition to finite current occurs when the mobil- 
ity edge \i c coincides with the vortex chemical potential 
/i. This is shown in more detail in Fig. 22, where the 
current is given as a function of the tilt for quarter filling 
n = 0.25. Apart from finite-size rounding of the tran- 
sition, the onset of the current is linear in h, consistent 
with the theoretical prediction in Eq. (80). 



VIII. SUMMARY 



predictions to numerics on finite systems which we di- 
agonalized exactly. We found good agreement with the 
RSRG results, although the numerics indicates that the 
universal regime is only accessible for very large system 
sizes. 

Our results seem to indicate that scaling pictures for- 
merly held regarding the case of uniform pins at a random 
separation are not valid near critical tilting, particularly, 
we find B± cx (H± — Hj_) in 1 + 1 dimensions, but with 
C = 1 rather than the expected £ = 1/2 (see Sees. I and 
VI D). We hope that our results will encourage further 
experimental and theoretical research on vortex pinning, 
both to verify our predictions, and to explore how strong 
randomness physics may appear in situations not touched 
upon here. 
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In this paper we studied the physics of interacting vor- 
tices in a 2-dimensional type-II superconductor with par- 
allel random columnar defects in the plane. Our analysis 
is based on mapping the system onto a 1+1-dimensional 
ensemble of hard-core bosons with on-site or hopping dis- 
order. The equivalence of hardcore bosons and fcrmions 
in one spatial dimension allowed us to focus on an ef- 
fectively noninteracting, but disordered, many-fermion 
system for the special case of Luttinger liquid parame- 
ter g = 1. As a qualitative benchmark, we considered 
the exactly solvable Lloyd model with on-site disorder 
where we obtained analytic predictions for the critical 
tilt and the transverse magnetic flux above the critical 
tilt. We then studied extensively the case of general ran- 
dom hopping, using the real-space renormalization group 
(RSRG) technique, which we generalized to the case of 
finite tilt. The RSRG allowed us to completely describe 
the localized phase and to extract physical quantities like 
the density of states or the mobility edge as a function of 
tilt. With the help of a non-Hcrmitian spectral formula 
we were able to extend the RSRG results to the delo- 
calized states, and to calculate the onset of transverse 
magnetic flux close to the breakdown of the transverse 
Meissner effect, which is the most relevant experimen- 
tal observable. Finally, we have compared our analytic 



APPENDIX A: A GENERAL FORMULA FOR 
THE TRANSVERSE VORTEX DENSITY NEAR 
THE MOBILITY EDGE 

If we increase the magnetic field above the critical field, 
some vortices are delocalized, and a transverse vortex 
density appears. The transverse vortex density (imagi- 
nary current) for a single vortex state is given by: 



J n = yRo— — - 
L oh 



(Al) 



Following Ref. 8,9, we write the eigenvalues of the ran- 
dom hopping problem as 



e n (h) = E(h + ik n ). 



(A2) 



The function E(h + ik) is analytic and therefore obeys 
the Cauchy-Riemann equations. In particular: 



„ dE dE 



(A3) 



From Eq. (A3) we can get a general result for the total 
current: 
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Kmax 

J tota; = 2 i^Re^i^2 J glm|i = i(Ini£ femax -Im^ min ) = i(Im^-Im^ m ) (A4) 



From Eq. (A4) we can immediately conclude, for in- 
stance, that closed bubbles of tilted states don't con- 
tribute to the total current. The factor of two accounts 
for the two branches of the spectrum. 

APPENDIX B: THE RSRG AND CRITICAL TILT 
FROM THE DISCRETE NON-HERMITIAN 
SPECTRAL FORMULA 

The discrete non-Hermitian formula allows an alter- 
native formulation of the real-space RG, and thus gives 



The real-space RG is obtained as follows. Assuming 
strong disorder, as we did in Sec. IV, we know that the 
eigenvalues ±ei are associated with the strongest bond 
of the chain, m: 

ei w w ni (B3) 

Let us choose an energy \E\ <C ei, and assume that it 
belongs to a localized state. We can now write: 

N/2 N 

£ In \E 2 - e 2 | = N\h\ + ^ In ( Wi ) In e 2 . (B4) 
i=2 i=l 

However, we can rearrange terms on the right hand side 
as follows: 

lnw ni -i + lnw ni + lnio ni+ i - lnef 

„ 1 / ^-ii^m + A _ i eff ( B5 ) 
~ m V w ni J — mW n 1 -l,n 1 +2' 

where w e J^_ x „ i+2 is just the effective hopping according 
to the real space RG as found in Eq. (39). We now 
have the renormalized chain with the two sites m, n\ + 1 
removed, and we can write the spectral formula as: 

N/2 N-2 

^2\n\E 2 -ej\ =N\h\ + ^ \n(w t )-\ne 2 . (B6) 

i=2 i=l 

where the w), are the hoppings of the renormalized chain. 
This process can now be repeated - starting with the 



I 

further support to the results obtained in Sec. VI. We 
will use the the real part of the non-Hermitian spectral 
formula, Eq. (70): 



N N 

^ - In ((RcE - e t f + (ImE) 2 ) - N\h\ + ^ In (w t ) . 
i=l i=\ 

(Bl) 

Because of particle-hole symmetry, we can also write: 



(B2) 

I 

next eigenvalue, 62 of the largest renormalized bond, 
max{wi} = w n2 , the bond ri2 can be renormalized in the 
same way as written above. Thus we recover the real- 
space RG - we can eliminate the high energy states and 
write the spectral formula in terms of the new couplings, 
and the reduced length. Note that the tilt part remains 
constant in the process. This invariance corresponds to 
the second condition in Eq. (39), which can be written 
as: 

* i 

We can repeat this process until we reach energy eigen- 
values of the same order as E. 

We can also re-derive the critical tilt result, Eq. (60). 
We assume that when close to the critical tilt, or to the 
mobility edge, the eigenvalue E, which we assume is e m < 
\E\ < e m+ i, is as far from e m , e m +i as it can be [such 
that the LHS of Eq. (B6) is maximized], and therefore: 

l]Tln||£|- ei |«l. (B8) 

where i <~ m indicates summing over i's in the vicinity 
of m and m + 1. Eq. (B8) is justified by the above 
assumption, and our understanding that in the random 
singlet phase, when renormalizing a bond with energy 
E, there is negligible probability to find another bond of 
similar strength next to it. With Eq. (B8), we can now 



N/2 N 

~ (In ((Rc£ - e t ) 2 + (ImE) 2 ) + In ((ReE + e t ) 2 + (ImE) 2 )) = N\h\ + ]T In ( Wl ) 



i=i 



i=i 
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write the spectral formula assuming we can carry out the 
RG procedure above until we arrive at the energy scale 
\E\. Also, since we ignore eigenvalues of the order of 
\E\, all remaining eigenvalues are much smaller than \E\. 
Therefore we can write: 

N/2 N 

^ln|£ 2 | KN\h\+^2]n(wi), (B9) 

i=l »=1 

where N is the length of the chain renormalized to energy 
scale \E\. This, in turn, can be written as: 

VlnM = N(T E +T a ) ^Nh c (BIO) 

where we used Eq. (46) to express the average of the log- 
arithmic couplings in terms of Te- We have now made 



the replacement \h\ — > h c , the critical tilt. Before, we as- 
sumed that the tilt is as close to critical as possible, and 
therefore the LHS of Eq. (BIO) is the maximum possi- 
ble value for N\h\, without letting E have an imaginary 
part. To get \h\ < h c we only need to allow E to come 
arbitrarily close to e m . Thus we conclude: 



N (Te + r ) = fp 

n rv + ry 



where we used Eq. (49) for the effective density of 
the renormalized chain. This concludes an independent 
demonstration of Eq. (60). 
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